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We present a coordinate-independent method for extracting mass (Ma) and angular momentum 
( Ja) of a black hole in numerical simulations. This method, based on the isolated horizon framework, 
is applicable both at late times when the black hole has reached equilibrium, and at early times 
when the black holes are widely separated. Assuming that the spatial hypersurfaces used in a 
given numerical simulation are such that apparent horizons exist and have been located on these 
hypersurfaces, we show how Ja and Ma can be determined in terms of only those quantities which 
are intrinsic to the apparent horizon. We also present a numerical method for finding the rotational 
symmetry vector field (required to calculate Ja) on the horizon. 



PACS numbers: 04.25.Dm 



I. INTRODUCTION 

Given a numerical simulation of a spacetime containing 
a black hole, one is faced with two immediate questions: 
Where is the black hole, and what are the values of its 
parameters, such as mass and angular momentum? 

In analytical considerations, a black hole is often de- 
fined via the event horizon, i.e. the future boundary of 
the causal past of future null infinity. While this notion 
is mathematically elegant and has led to powerful results 
in black hole physics, generally it is not directly useful in 
numerical evolutions of black hole spacetimes. The rea- 
son is the teleological nature of the event horizon: it can 
be constructed only after we have knowledge of the entire 
spacetime and spacetime is the end product of these sim- 
ulations. Furthermore, due to practical limitations, gen- 
erally one cannot evolve all the way to future null infin- 
ity. For such reasons, in numerical evolutions, it is more 
common to use apparent horizons to characterize a black 
hole. Apparent horizons are closed two-surfaces on a spa- 
tial slice of the spacetime and are, therefore, well suited 
to numerical simulations that evolve data from one spa- 
tial slice to another. While it is possible that there exist 
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no apparent horizons on the spatial slices used in a given 
numerical simulation of a black hole spacetime 1 , in prac- 
tice, this issue does not arise for most numerical black 
hole simulations that we are aware of. We shall therefore 
disregard this issue for the rest of this paper and only deal 
with partial Cauchy surfaces (in black hole spacetimes) 
for which apparent horizons do exist and have been lo- 
cated 0, II H II i, H H- A natural question we want 
to ask is: can one unambiguously associate black hole 
parameters to apparent horizons? 

One way to attribute a mass and an angular momen- 
tum to a black hole is to calculate the corresponding 
ADM quantities at infinity. The main difficulty is that 
the ADM mass and angular momentum refer to the whole 
spacetime. In a dynamical situation, such a spacetime 
will contain gravitational radiation and it is not clear 
how much of the mass or angular momentum should be 
attributed to the black hole itself and, if there is more 
than one black hole, to each individual black hole. 

It is desirable to have a framework that combines the 
properties of apparent horizons with the powerful tools 
available at infinity. In the regime when the black hole 
is isolated in an otherwise dynamical spacetime, such a 
framework now exists in the form of isolated horizons 
El EH E3l • In numerical simulations of, say, black 
hole collisions, the black hole would be isolated at early 
times when the black holes are well separated, or at late 
times when the final black hole has settled down, but ra- 
diation is still present in the spacetime. Isolated horizons 



For example, it has been shown in |l| that there exist Cauchy sur- 
faces containing no apparent horizons even in the Schwarzschild 
spacetime 
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provide a way to identify a black hole quasi-locally, and 
allow for the calculation of mass and angular momen- 
tum. It has been shown recently that the formulae for 
angular momentum and mass arising from the isolated 
horizon formalism are valid even in dynamical situations 
[l4j . Therefore, the numerical methods presented in this 
article are relevant even for fully dynamical black holes. 

The aim of this paper is to introduce the isolated hori- 
zon framework in a way that is directly useful in nu- 
merical relativity. In particular, we show how to find 
isolated horizons numerically, and how to implement the 
isolated horizon formulae for Ja and A key result 

of this paper is that, to find the angular momentum of 
an isolated horizon, one can just use the ADM formula 
for angular momentum but now applied at the apparent 
horizon. This is not an assumption, but a rigorous re- 
sult obtained by calculating the Hamiltonian generating 
diffeomorphisms which reduce to rotational symmetries 
on the isolated horizon. This is completely analogous 
to what is done at infinity to obtain the ADM formulae 
for mass and angular momentum for asymptotically flat 
spacetimes. Indeed, the Hamiltonian analysis of isolated 
horizons is an extension of the ADM formalism to the 
case where the region of spacetime under consideration 
has an inner boundary in the form of an isolated hori- 
zon. The isolated horizon results for Ja and Ma are also 
convenient for practical reasons because their expressions 
only involve data defined on the apparent horizon. An 
important ingredient in the formula for Ja is an axial 
symmetry vector on the apparent horizon. Therefore we 
also present and implement a numerical method for lo- 
cating Killing vectors on the horizon. 

This paper is organized as follows. In section [n] we 
first review the notion of isolated horizons and explain 
how they can be located in numerical simulations of black 
hole spacetimes. We also explain their relation to appar- 
ent horizons and give the formulae for their mass and 
angular momentum. These formulae require that we find 
a Killing vector on the apparent horizon. A method for- 
doing this is described in section lTTTl Scction llVl describes 
the practical numerical implementation of these formu- 
lae, and we verify that our method gives the expected 
results in the simple case of a boosted Kerr black hole. 
Section |V| compares the isolated horizon formulae with 
other approaches used to compute mass and angular mo- 
mentum, and finally in section IvTI we give some possible 
applications and future directions. 



II. ISOLATED HORIZONS 

Let us first fix our notation. The spacetime metric is 
taken to have signature (— , +, +, +), and we use geomet- 
rical units where G and c are equal to unity. We will 
usually use the abstract index notation, but occasion- 
ally, especially for differential forms, the index free no- 
tation will also be used. The Riemann tensor R a bc d will 
be defined by the equation 2V[ a Vb]0! c = R a bc d &d where 



a a is an arbitrary co-vector. We work on a spacetime 
that is foliated by a spatial three manifold £ which is 
a partial Cauchy surface. The spacetime is thus of the 
form Ext. The three-metric and extrinsic curvature 
on a spatial slice £ are respectively denoted by jab and 
K ao := — 7a c 7b d V c Td where T a is the unit timelike nor- 
mal to S. The two- metric on the apparent horizon is 
called q a b. All manifolds and fields will be taken to be 
smooth. Finally, though it is quite easy to include mat- 
ter fields, we will be interested in vacuum spacetimes 
only. This paper is reasonably self contained, however, 
for precise and detailed definitions and proofs we refer the 
reader to 0, 0, 0, ^| . In this paper we shall focus on 
the basic physical ideas underlying the isolated horizon 
framework and its application to numerical relativity. 

A. Motivation, definition and basic results 

For completeness and to fix notation, let us start by 
reviewing the definition of apparent horizons. Let T a 
be the unit timelike vector field orthogonal to a partial 
Cauchy surface S. Given a closed two-surface S C S, we 
have the unique unit outward-pointing spacelike normal 
R a which is tangent to E. We are assuming here that 
one can define what is meant by the inside and outside 
of S. This can be done, for instance, by imposing cer- 
tain global conditions on the spacetime and the partial 
Cauchy surfaces under consideration. These conditions 
include: asymptotic flatness, strong future asymptotic 
predictability of the spacetime with the Cauchy surfaces 
under consideration and certain topological restrictions 
on E. Since the focus of this paper is not on such global 
topological properties of spacetimes, we shall not spell 
out these assumptions precisely and instead refer the 
reader to [l5| . It is worth mentioning that the inside and 
outside S can also be defined quasi-locally without using 
global assumptions such as the ones mentioned above. 
This has been done by Hayward 0] using the notion of 
trapping horizons which we discuss in a later section. 

Let q a b be the induced Riemannian two-metric on S 
and e a b the area two- form on S constructed from q a i,. 
We can construct a convenient basis for performing cal- 
culations at points of S in a natural way (see figure 
First, define the outgoing and ingoing null vectors 

r := -Ur a + R a ) and n a := ^=(T a - R a ) . (1) 

It is worth noting that any spacelike two-surface S deter- 
mines uniquely, up to rescalings, two null vectors orthog- 
onal to S. Any other choice of t and n will differ from the 
one made in equation Q only by possible rescalings. We 
tie together the scalings oi£ and n by requiring £-n = — 1. 

Next, given two arbitrarily chosen spacelike orthonor- 
mal vectors e\ and e 2 tangent to S, construct a complex 
null vector 

m := —7={ei + ie 2 ). (2) 
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from the definition of an apparent horizon, the outward 
normal I is always expansion free. 

We will now argue that the isolated horizon frame- 
work is ideally suited to describe apparent horizons 
in the regime when the world tube H is null. For 
our purposes, the straightforward definition of a non- 
expanding horizon given below shall suffice. To carry 
out the Hamiltonian analysis in order to define mass and 
angular momentum, we actually need to impose further 
conditions on non-expanding horizons, which we shall 
briefly describe towards the end of this section. The 
formulae for mass and angular momentum make sense 
even on non-expanding horizons. 



FIG. 1: The figure shows an apparent horizon S embedded 
in a spatial slice E. T a is the unit timelike normal to E and 
R a is the outward pointing unit spatial normal to S in E; l a 
and n a are the outgoing and ingoing null vectors, respectively. 
The vector m a (not shown in the figure) is tangent to S. H 
is the world tube of apparent horizons. 



It satisfies the relations m ■ m = 0, m ■fn = 1, £ ■ m = 0, 
and n ■ m = 0. Since I and n satisfy £ ■ n = —1, we 
see that (£, n, m,m) form a null tetrad at 5*. This is, 
of course, only one possible choice of null tetrad, and 
we must ensure that physical results are independent of 
this choice. The expansions of £ a and n a are defined as 
0(i) :— q ab X7 a £b and 6r n ) := g o6 V a n; ) , respectively. Note 
that in order to find the expansions, we only need deriva- 
tives of £ and n along S, and there is no need to extend 
the null tetrad into the full spacetime. However, if in 
some numerical computations it is necessary to extend 
the null tetrad smoothly into the full spacetime; all cal- 
culations will be insensitive to this extension. The surface 
S is said to be an apparent horizon if it is the outermost 
outer-marginally-trapped-surface, i.e. it is the outermost 
surface on £ with 0(t) — and 0i n \ < 0. 

Consider now the world tube of apparent horizons H 
constructed by stacking together the apparent horizons 
on different spatial slices. As we shall show later in sec- 
tion ^TOJ smooth segments of this world tube are gener- 
ically spacelike; it is null when no matter or radiation is 
falling into the black hole. At late times, one expects the 
black hole to reach equilibrium when radiation and mat- 
ter are no longer crossing the horizon. In this regime, the 
world tube H will be a null surface, and the two-metric 
q a b on the apparent horizon S may now be viewed as a 
degenerate three-metric on the null surface H . Further- 
more, from experience with numerical simulations and 
also from general topological censorship results (see e.g. 
H3)> we know that at late times, the apparent horizons 
must have spherical topology. Therefore, at late times, 
the topology of H is S 2 x R. Finally, in this regime, 
the outward-normal I constructed in equation P is a 
null normal to the world tube H, and most importantly, 



Definition: A three dimensional sub-manifold A of a 
space-time (M,g a b) is said to be a non-expanding horizon 
(NEH) if it satisfies the following conditions: 

(i) A is topologically S 2 xl and null; 

(ii) The expansion 9 m := q ab V ' a £b of I vanishes on A, 

where I is any null normal to A; 

(iii) All equations of motion hold at A and, if any mat- 
ter field are present with T a b as the stress energy 
tensor, then we require —T£i b to be future directed 
and causal for any future directed null normal £ b . 

Note that if condition (ii) holds for one null normal I, 
then it holds for all. We will only consider those null nor- 
mals which are nowhere vanishing and future directed. 
We are therefore allowed to rescale £ by any positive- 
definite function. The energy condition in (iii) is implied 
e.g. by the null energy condition which is commonly as- 
sumed; for the remainder of this paper, we shall usually 
consider only vacuum spacetimes. 

Comparing the properties of the world tube H de- 
scribed earlier with conditions (i) and (ii) in the defi- 
nition, we see that the NEH is precisely what we need to 
model the physical situation at hand; when the black hole 
is approximately isolated, the world tube H represents a 
non-expanding-horizon A. The motivation behind the 
conditions in the definition are thus rather straightfor- 
ward from the perspective of apparent horizons. 

We have motivated the above definition using the world 
tube of apparent horizons which in general clearly de- 
pends on how we choose to foliate the spacetime with 
spatial surfaces. It is then natural to ask whether an 
NEH depends on the choice of this foliation. To avoid 
unnecessary technicalities, let us consider a possibly dy- 
namical spacetime containing one isolated black hole. As 
mentioned in the introduction, we can choose spatial hy- 
persurfaces on which there are no apparent horizons so 
that the world tube of apparent horizons does not exist 
for this foliation. Let us ignore this possibility and fo- 
cus on the case when each spatial slice has exactly one 
marginally trapped surface and hence exactly one appar- 
ent horizon; the generalization to multi black hole space- 
times is straightforward. It is not difficult to verify that 
the world tube of apparent horizons is independent of 
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the spacelike hypersurfaces used to construct it. In other 
words, the null surface (the NEH) obtained by stacking 
together all the apparent horizons at different times is 
independent of how we choose to foliate the spacetime 
by spacelike hypersurfaces; this is not true in the regime 
when the horizon is non-isolated and the world tube is 
spacelike. An NEH, even when viewed as the world tube 
of apparent horizons, is therefore an invariant notion in 
the four dimensional spacetime. Every cross section of an 
NEH is potentially an apparent horizon if it arises as the 
intersection of the NEH with a spatial slice. Choosing a 
different spatial slice corresponds to obtaining a different 
cross section of the NEH. As we shall see later, for cal- 
culating physical quantities such as angular momentum 
and mass, it does not matter which cross section of the 
NEH we choose to work with. 

Even though cross sections of an NEH are very closely 
related to apparent horizons, there are two differences: 
(i) Apparent horizons are required to be the outermost 
surfaces on a spatial slice with the afore-mentioned 
properties. This may not be true in general for cross 
sections of A; (ii) Since they are trapped surfaces, 
apparent horizons also satisfy the condition Oua < 0. 
Though this will most likely be true in actual numerical 
simulations, it turns out that this condition is not 
required to study the mechanics of isolated horizons. 
In fact, there exist exact solutions representing black 
holes which are isolated horizons but do not satisfy 
0i n ) < 0, e.g., the distorted black holes studied by 
Geroch and Hartle ^^|. In these solutions the inte- 
gral of 6>(„) over a cross section of the horizon is still 
negative even though 0r n ) is not necessarily negative 
everywhere. In the remainder of this paper, we shall 
ignore these caveats and the phrases 'apparent hori- 
zon' and 'cross-section of A' will be used interchangeably. 

Consequences of the definition: The simple definition of 
a NEH leads to a surprisingly large number of conse- 
quences. The important results we shall need later are 
the existence of a one-form lu (1 whose definition is given 
in equation and the fact that the component ^2 of 
the Weyl tensor is a physical, i.e. gauge invariant, prop- 
erty of an isolated horizon. Keeping these two results in 
mind, the rest of this sub-section may be skipped on a 
first reading. 

While stating the results, it is convenient to use a null- 
tetrad (£, n, m,m) which is adapted to the horizon] this 
means that £ is a (future directed) null normal to A. An 
example of such a null-tetrad is the one constructed in 
equations Q and @, but there is, of course, an infinity of 
such tetrads. Physical quantities will be independent of 
which null-tetrad we choose. These consequences do not 
require the spacetime to be vacuum; they only assume the 
energy condition mentioned after the definition above. 

• Area of apparent horizons is constant: Any null 
normal £ is always geodetic and is also a symme- 
try of the degenerate intrinsic metric q a b on A: 
£ eqab — (the symbol '=' means that the equality 



holds only at points of A). Furthermore, the area 
of any cross-section S defined as A a ■— § s d 2 V is 
the same where d 2 V is the natural area measure on 
S constructed from the two- metric q a b- In partic- 
ular, the area of the apparent horizons is constant 
in time. 

• Definition of uj a : It will be shown in section 111 HI 
that the shear of £ defined as :— m a m b V a £b is 
zero. The twist of £ is trivially zero because I is 
normal to a smooth surface, and the expansion of 
£ is zero by definition. Since the null normal I is 
expansion, twist and shear free, there must exist a 
one-form uj a intrinsic to A such that 

t a V a £ b 4 t a u>J b (3) 

for any vector t a tangent to A. If the null normal I 
is rescaled I — > £ = f£, the one-form u> a transforms 
as 

uj — * uj = uj + d (In/) (4) 
where 'd' is the exterior derivative on A. 

• Gauge invariance of ^2' If we choose a null tetrad 
adapted to A such as in equations and @ , then 
we get the following restrictions on the Weyl tensor 
Cabcd at A [HI: 

*o = *i = (5) 

and 

duj = 21m [# 2 ] e, (6) 

where 

*o = C abcd £ a m b £ c m d (7) 
*i - C abcd £ a m b £ c n d (8) 
*2 = C abcd £ a m b m c n d (9) 

and e is the area two-form on the apparent horizon. 
This can be used to prove the important result that 
^2 is gauge invariant. The gauge freedom we are 
concerned with here is the choice of null tetrads 
adapted to A. The allowed gauge transformations 
relating different choices are null rotations about b. 

£ -» £ 

m — > m + c£ (10) 
n — > n + cm + cm + cc£ 

and spin-boost transformations: 

£ -> A£ 

n -> A^n (11) 

m — > e 2ie m 

Here c, A and 9 are arbitrary smooth functions on 
A. It turns out that ^2 is always invariant under 
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spin-boost transformations; under null rotations, it 
transforms as (see e.g. the appendix of |19|) 

^ 2 * 2 + 2c*i +c 2 * - (12) 

Since ^>o and \&i vanish at A, we see that 1 $>2 is 
in fact gauge invariant at the horizon: it does not 
depend on the choice of null tetrad as long as £ is 
one of the null generators of A. This property is 
important because it tells us that the value of \&2 
at the horizon does not depend on how we choose 
to foliate our spacetime. A different spatial slice 
£ will lead to a different £, n and m. The two 
null tetrads will be related by a combination of the 
following transformations: a null-rotation about £, 
a spin-boost transformation, or a multiplication of 
m by a phase. Whichever null-tetrad we use to 
calculate ^2, we will get the same result. As we 
shall see, it is the imaginary part of \&2 which is 
physically interesting for our purposes. 

Further Definitions: The notion of non-expanding hori- 
zons describes the late time behavior of apparent hori- 
zons. In order to define the mass Ma and angular mo- 
mentum Ja of A, one needs to go beyond this defini- 
tion and introduce additional structures on the horizon. 
This is done via the definition of a weakly isolated hori- 
zon 0, 0| . The Hamiltonian analysis which leads to 
the definitions of mass and angular momentum requires 
this extra structure. Fortunately, it turns out that the 
formulae for Ma and Ja do not depend on this extra 
structure and hold true for non-expanding horizons. We 
could therefore omit these definitions entirely and sim- 
ply state the results of the calculation. For the sake of 
completeness, we shall give the basic idea behind weakly 
isolated horizons. The rest of this sub-section may be 
skipped without loss of continuity. 

In an NEH, the intrinsic metric q a b is time indepen- 
dent since C iq a b — 0. There is no restriction on the time 
derivatives of the extrinsic curvature of A or the intrinsic 
connection on A (by 'time derivative' we mean derivative 
along £). Since A is a null surface, there is no natural no- 
tion of extrinsic curvature. The closest thing to extrinsic 
curvature is the tensor K a defined via t a K a b = t a \7 a £ b 
for any t a tangent to A. This tensor is known in the 
mathematics literature as the Weingarten map, and as 
its definition shows, it is V a £ fc with the covariant index 
pulled back to A. From equation @ we immediately ob- 
tain K a = uj a £ b . This implies that requiring K a to be 
time independent is equivalent to requiring Ctu a = 0. 
Unfortunately, due of the transformation property of uj a 
(see equation Q), it is clear that this equation is not 
meaningful if all rescalings of £ are allowed. Note that 
if we restrict ourselves to rescalings that are constant on 
the horizon, then to a is invariant. We thus need to restrict 
ourselves to an equivalence class of null normals [£] , the 
members of which are related to each other by a con- 
stant, positive non-zero rescaling. We can then associate 
a unique u a with [£] . The equation C £0J a = is now 



perfectly meaningful if I is a member of [£]. A weakly 
isolated horizon is then defined as an NEH A equipped 
with such a preferred equivalence class of null normals [£] . 
Furthermore, the Lo a associated with [£] must be time in- 
dependent: C iuj a = 0. Given an NEH, we can always 
find such a preferred equivalence class (but this equiva- 
lence class is not unique) . Thus every NEH can be turned 
into a weakly isolated horizon. The Hamiltonian analy- 
sis of weakly isolated horizons leads to formulae for mass 
and angular momentum, and as we shall see in section 
III CI these formulae are insensitive to arbitrary rescalings 
of £ and thus they make sense even for non-expanding 
horizons. 

We can introduce an even stronger definition. A 
weakly isolated horizon requires that co a be time inde- 
pendent. It can be seen from equation © that Lu a is 
also a component of the intrinsic connection D a induced 
on A by the four dimensional connection V a compatible 
with the four-metric. However, u a is just one compo- 
nent of D a . In an isolated horizon, we require that all 
components of D a are time independent. It turns out 
that generically, this condition gives us a preferred unique 
equivalence class [£] |l3| . 

B. Finding non-expanding horizons 

Our strategy to find non-expanding horizons is to lo- 
cate apparent horizons on each spatial slice and then to 
check whether the world tube H obtained by stacking 
these horizons together is a NEH. What remains to be 
checked is whether the tube is a null surface. By defini- 
tion, a surface is null if the metric h a b induced on this 
surface has a degenerate direction, i.e. if there exists a 
vector X a tangent to H such that h a bX b = 0; therefore, 
one possible method to check whether H is isolated is 
to construct the induced metric h a b on H and see if it 
has a zero eigenvalue. To construct h a b numerically, we 
have to know the two-metric q a b on at least two differ- 
ent time slices. Furthermore, in a numerical simulation, 
H will never be exactly isolated because of numerical 
errors, and it is not clear how this method can quan- 
tify how close the horizon is to being exactly isolated. 
Fortunately, there is a much simpler method which only 
requires the two-metric on a single time slice and can 
also provide a quantitative measure of how close H is to 
being perfectly isolated. This method is based on the 
shear o b of £, which is the symmetric trace-free part 
of the projection of \7 a £b onto the apparent horizon S. 
The tensor <r ab has two independent components, and is 
conveniently written in terms of a single complex num- 
ber :— m a m b W a £b, where m is defined in equation 
© ■ To calculate onn conveniently, we simply decompose 
£ using equation Q: 

= m m V a £b 

= -^m a m b \7 a T b + -^m a m b \7 a R b . (13) 
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In the above equation, we have extended the notation 
'=' to also mean that the equality holds only at points 
of H. The first term is just a component of the ex- 
trinsic curvature K a f )1 while the second term can be 
calculated on the spatial slice by calculating the con- 
nection associated with the three-metric 7 a &. We shall 
now prove the following important result concerning a^y 

The world tube of apparent horizons is a NEH if and 
only if am = 0. 

Readers not interested in the proof of this statement 
can skip forward to equation (|2U|) . 

To prove this statement, we need to consider the gen- 
eral case when H is not null. This has been studied in 
great detail by Hayward |l6j |. In Hayward's terminol- 
ogy, the surface H is essentially a future outer trapping 
horizon. This means that H is foliated by a family of 
marginally trapped surfaces (which in our case are the ap- 
parent horizons) satisfying the relations 9 m = 0, 9r n \ < 
and C n 9m < 0. These are physically very reasonable 
conditions, and all black holes found in simulations are 
expected to satisfy them. Incidentally, note that the con- 
dition C nfiu) — can be used to define, in a quasi local 
manner, what we mean by the inside and outside of the 
horizon. 

The proof of this statement, adapted from is then 
quite simple: let z a be a vector tangent to H and or- 
thogonal to the foliation whose leaves are the apparent 
horizons. Such a vector field may be considered to define 
time evolution at the horizon. It is easy to see that, up to 
a rescaling, z a can be expressed as a linear combination 
of £ a and n a 



l a 



an 



(14) 



where a is a smooth function on H. The rescaling free- 
dom in z will be inconsequential for our purposes. The 
surface H is null, spacelike, or timelike if and only if a is 
zero, positive, or negative respectively. We will now show 
that a > 0. From the definition of apparent horizons we 
know that the expansion 9m vanishes everywhere on the 
horizon, therefore C z 9m = 0. This in turn gives 



a 



(15) 



Even though 9m and 9i n ) are so far defined only on iJ, 



in equation (|15|) (and also in the very definition of a trap- 
ping horizon) we are taking the derivatives of these quan- 
tities along £ and n which are not necessarily tangent to 
H. To make sense of this equation we need to extend 
£ and n in a neighborhood of the surface S. This can 
easily be done by using the unique geodesies determined 
by these vectors. 

The Raychaudhuri equation for I then leads to 



C 



(16) 



where $oo = \R a b£ a £ h ■ If we assume that the spacetime 
is vacuum, then $oo = 0, and therefore C i9m — ~\ a U) | 2 > 
which along with equation l|15|) gives 
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c 



(17) 



This immediately implies that a = (which is equivalent 
to H being null) if and only if am = 0. This is what we 
wanted to show. More generally, since C n 9m < 0, this 
shows that a > 0, which means that A is spacelike when 
the shear is non-zero. 

As a side remark we also show that H is null if and 
only if the area element on the apparent horizons e a b is 
preserved in time. To show this we need the equations 

C-itab = 9{i)e a b = and £„e a fc — 9( n )^ab ■ (18) 
It then follows that 



•£ z^ab — £ ie a 



i^ab 



-a.9. 



(n)£ab 



(19) 



therefore a = if and only if C z e a b — 0, which is what we 
wanted to prove. This implies that if a = 0, then the area 
of cross-sections of A is constant. However, the converse 
is not necessarily true, because e a b could be changing in 
such a way that its integral is constant; the area can be 
constant globally without being constant locally. Thus, 
in principle, we can have situations in which the area is 
constant without the horizon being isolated. However, 
constancy of area is still a very useful first check to see 
when the horizon reaches equilibrium. Finally, as a side 
remark we note that since a > and 6i n ) < 0, we get 
£ z Aa > 0, which is the area increase law. 

In this paper we are interested in the case when am 
vanishes (up to numerical errors). In order for the hori- 
zon to be isolated we should have 



2 d 2 V 







(20) 



where d V is the natural area measure on S constructed 
from q a b. The quantity s is dimensionless since am has 
dimensions of inverse length. For the horizon to be nu- 
merically isolated, we require that s converge to zero ap- 
propriately when the numerical resolution is increased 2 . 



2 We want to point out that as it stands, the quantity s can not 
be used as a general measure of how isolated a given horizon is, 
since it is not gauge invariant. If s is not exactly zero to begin 
with, a simple rescaling of I changes cr^ and thus s. However, 
if we are only concerned with given apparent horizons embedded 
in given spatial slices, and if we agree to use equation Q for 
defining I thereby removing the boost freedom, then the horizon 
will be close to being isolated if the condition s < 1 is satis- 
fied. While this is not a satisfactory solution for identifying the 
small parameter describing a horizon which is 'almost isolated', 
it might nevertheless be useful in practice. The com plet e solu- 
tion to this problem will require using results from 14] where 
the world tube of apparent horizons is analyzed for the spacelike 
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C. Mass and angular momentum 

Now let us define angular momentum. The details and 
proofs of the results below can be found in . In order 
to define the angular momentum, we need to assume that 
S is axisymmetric, i.e. there is a vector field ip a tangent 
to S such that it preserves q ab 3 

£ v q ab = 0. (21) 

The next section describes a simple method for finding 
(p a . Note that ip a is not a Killing vector of the full space- 
time metric. It is completely intrinsic to A, and in fact 
we only need data on the apparent horizon to find it. 

Given ip a and ui a (defined in equation the isolated 
horizon angular momentum Ja is defined as the surface 
term at A in the Hamiltonian which generates diffeomor- 
phisms along a rotational vector field which is equal to 
tp a at the horizon [13J. The expression for Ja is: 



Ja 



1 



(uj a ^ a )d 2 V 



= -— & /Im[tf 2 ] d 2 V, 
•J s 



(22) 



where S is the apparent horizon and the function / is 
related to ip a by d a f — £ ba tp b ■ In the second equality, we 
have used equation © and an integration by parts. Since 
^2 is gauge invariant, so is the angular momentum, and 
in particular, it does not depend on the scaling of I and 
it thus makes sense even on a NEH. If, in a neighborhood 
of A, there were a spacetime rotational Killing vector cf> a 
that approached tp a at A, then the above formula for Ja 
would be equal to the Komar integral calculated for cf> a 
[l2|. However, the formula in equation (|22|l is more gen- 
eral. As a practical matter, even if there were a Killing 
vector in the neighborhood of A, it is easier to use 1221) 
rather than the Komar integral, since the Komar inte- 
gral requires finding the Killing vector in the full four- 
dimensional spacetime. It is also worth mentioning that, 
if the vector field ip a is not a symmetry of A but is an 
arbitrary vector field tangent to S, then Ja is still the 
surface term at A in the Hamiltonian generating diffeo- 
morphisms along vector fields which are equal to ip a at 
the horizon. But in this case, there would be no reason 
to identify Ja with the angular momentum, since con- 
served quantities such as mass and angular momentum 



3 This is actually an oversimplification: in addition to preserving 
q a b, the vector field ip a must also preserve the additional struc- 
ture introduced to define weakly isolated horizons (discussed 
briefly at the end of section III Al . namely the preferred equiva- 
lence class [£] and the one-form u) a associated with [£] : C ipU! a = 
and C v £ = cl where c is a positive-definite constant. However, 
these additional conditions are easy to verify in practice, and 
the non-trivial condition is equation 1211 . We shall not concern 
ourselves with these extra conditions in the remainder of this 
paper. 



are always associated with symmetries. However, the ex- 
istence of the axial symmetry is the least that must be 
true if the horizon is to be close to Kerr in any sense. 
Note that angular momentum is a coordinate indepen- 
dent quantity; even if we use corotating coordinates to 
describe the black hole, the black hole still has the same 
angular momentum. 

We now describe a form of equation (1221 which is much 
better suited for calculating Ja numerically. From equa- 
tion (|3J|, which is the defining equation for uj a , we get 



t a w a ± -n b t a w a e b 4 e»t a v a n b 



(23) 



where, as before, t a is any vector tangent to A. Assume 
that we have found the symmetry vector field <p a on the 
horizon (the method we use for finding tp a is described 
below). From equation (|22|l . we are eventually interested 
in calculating <p a u> a ', since ip a is tangent to A, setting 
t a = (p a we get 

<p a w a 4 tp a tv a n b 

± ^ a (T b + R b )V a (T b -R b ) 

— \ V a {T b V a T b -T b V a R b 
+R b V a T b - R b V a R b ) 



(24) 



ip a R h V a T b 4 -p a R b K ab 



where we have used equation Q along with the fact that 
A^ a and R a are orthonormal. In the last step, the defini- 
tion of extrinsic curvature K ab = —J a c J b d VcTd has been 
used where ^ a b = 9ab + T a T b is the three-metric on the 
spatial slice E. We have thus reduced the calculation of 
(p a Lu a t° finding a single component of the extrinsic cur- 
vature. The integration of this scalar over the apparent 
horizon yields the angular momentum: 



Ja = ^j>{^ a R b K ab )d 2 V. 



(25) 



This is our final formula for the angular momentum. This 
formula is remarkably similar to the formula for the ADM 
angular momentum computed at spatial infinity: 



J 



ADM 



^£ {K ab - lab K)cp a A 2 S b 



2 Q b 



— K ab rd'S 



(26) 



The ^ a bK term does not contribute, because <f> a is tan- 
gent to Soo, which is the sphere at spatial infinity. Since 
the metric on is just the standard two-sphere met- 
ric, we have no difficulty in choosing a (f> a , and we can 
calculate Jadm about any axis. In contrast, since the 
metric on the apparent horizon S is distorted, finding (j) a 
is more complicated. Finally, as mentioned earlier, the 
similarity between equations l|25|) and (|26|l is not sur- 
prising, because both quantities are defined to be surface 
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terms in Hamiltonians generating diffeomorphisms along 
the appropriate rotational symmetry vector fields; Ja is 
the surface term at A while j\ DM is the surface term at 
infinity. 

Given Ja , the horizon mass Ma is given by IT^ 
Ma = ^yjRi + ^J 2 A (27) 

where Ra is the area radius of the horizon: i?A = 
(Aa/^tt) 1 / 2 . This formula depends on i?A and Ja in 
the same way as in the Kerr solution. However, this is 
a result of the calculation and not an assumption. The 
formulae for Ja and Ma are on the same footing as the 
formulae for the ADM mass and angular momentum, be- 
cause both are derived by Hamiltonian methods. Fur- 
thermore, under some physically reasonable assumptions 
on fields near future time- like infinity (i + ), one can show 
that Ma — Madm is equal to the energy radiated across 
future null-infinity if the isolated horizon extends all the 
way to i + . Thus, Ma is the mass left over after all the 
gravitational radiation has left the system. This lends 
further support for identifying Ma with the mass of the 
black hole. Further remarks regarding equation (|27[1 ap- 
pear towards the beginning of sec. 

To summarize: in order to calculate the mass and an- 
gular momentum of an isolated horizon, we need the fol- 
lowing ingredients: 

1. We must find the apparent horizon and check if the 
shear of the outward null normal vanishes within 
numerical errors. If it does, then the isolated hori- 
zon formulae are applicable. 

2. We need to find the quantity K a bR b on S and the 
symmetry vector ip a (assuming it exists). 

3. The angular momentum is then a simple integral 
over the apparent horizon given by equation (|25|l . 
and the mass is a purely algebraic function given 
by equation (|T7Jl . 

The first and last steps are rather straightforward if we 
know the location of the apparent horizon. The only non- 
trivial step is the calculation of ip a . In the next section, 
we show how it is numerically calculated using the two- 
metric on S. 



III. FINDING THE KILLING VECTOR 

First of all, we should point out that in some numeri- 
cal simulations (especially simulations with built-in axi- 
symmetry) the axial symmetry vector is already known. 
In that case, one can go ahead and find the angular mo- 
mentum using equation (|25() . However, we are also in- 
terested in the more general case, when there is an axial 
symmetry, but the coordinates used in the simulation are 
not adapted to it. In this section, we describe a general 



numerical method for finding ip a . Our method of find- 
ing Killing vectors on the apparent horizons is based on 
the Killing transport equation, which we now describe. 
This method does not depend on the fact that we are on 
an apparent horizon, and it is possible that this proce- 
dure could find Killing vectors efficiently in more general 
situations. We first describe the general method. 

Let £° be a Killing vector on (S,q a b), and define the 
two- form L a b — V a £fc. This is a two-form because of 
the Killing equation V( a £{,) = L( Q {,) = 0. It is then not 
difficult to prove the following (see e.g. [2(|) 

V a Vatb = V a L ab 

and v a V a L bc = Rcb/CdV a . (28) 

The reason for inserting an arbitrary vector v a will soon 
become clear. Instead of viewing these as equations for a 
Killing vector, let us instead think of them as equations 
for an arbitrary vector £ a (or a one- form £ a ) and an ar- 

(p) 

bitrary two- form L a b- If we start with a one- form Q[ 
and a two-form L^? at a point p on the manifold, then 
the above equations can be solved along any curve j(t) 
(with v a as its tangent) starting at p to give a unique one- 
form rj a and a unique two-form a a b at any other point on 
the curve. This procedure is analogous to parallel trans- 
port, but the differential equation used in the transport 
is not the geodesic equation, but instead equation 128|) 
above, and instead of transporting a vector, these equa- 
tions transport a one-form and a two-form. Viewed this 
way, equations l|28|l are often referred to as the Killing 
transport equations [21|. We are thus led to consider the 
vector space V p consisting of all pairs (ffi , ) for an 
arbitrary one-form and an arbitrary two-form L-v 
at a point p. For any curve j(t) which starts at p and 
ends at q, the equations in (J2SJ, being linear, give us a 

linear mapping between V p and V q . If (^ } ,L^) £ V p 
actually comes from a Killing vector and its derivative, 
then it will be mapped to (£i 9 ' , V?) ) G V q , which comes 
from the same Killing vector. 

If we consider closed curves starting and ending at the 
point p, then the Killing transport for a curve j(t) gives 
us a linear mapping M p {^) : V p — ► V p . A Killing vector 
corresponds to an eigenvector of M p {-f) with eigenvalue 
equal to unity for any closed curve 7. Note that an eigen- 
vector Z of M p (7) with unit eigenvalue need not neces- 
sarily correspond to a Killing vector; this will be true 
only if Z is an eigenvector of M p (7) with unit eigenvalue 
for every sufficiently smooth curve 7 beginning and end- 
ing at p. However, if there exists only one eigenvector 
with eigenvalue 1, it is the only possible candidate for a 
Killing vector and it is easy to check explicitly (by calcu- 
lating the appropriate Lie derivative of q a b) whether this 
candidate is indeed a Killing vector. If there exists no 
such eigenvector, then it tells us that q a b does not have 
any Killing vectors at p. 

In our case, S is a topological two-sphere which means 
that V p is a three dimensional vector space and, if we 
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choose a basis, M p (-f) can be represented as a 3 x 3 ma- 
trix (if S is an n dimensional manifold, then Mp(j) is 
a \n{n + 1) dimensional matrix). Finding the Killing 
vector at p then reduces to an eigenvalue problem for a 
3x3 matrix. For a constant curvature two-sphere, as 
in the Schwarzschild horizon, this matrix will just be the 
identity matrix for any point p. For an axially symmet- 
ric sphere, such as the one in a Kerr spacetime, there 
will be precisely one such eigenvector. Having found £ a 
and L a b at one point, we can again use equation l|28|l 
to find it everywhere on the sphere, using various other 
curves. Finally, the Killing vector is normalized by re- 
quiring its integral curves to have affine length 2tt (it 
can be shown that the integral curves must in fact be 
closed). Since we are only free to rescale the Killing vec- 
tor by an overall constant, we only have to perform the 
normalization on one integral curve. This normalization 
is valid for rotational Killing vectors. If we were dealing 
with, say, translational or stationary Killing vectors, the 
appropriate normalization condition would be to require 
the vector to have unit norm at infinity. 

To make this procedure concrete, let us write down the 
equations explicitly in spherical coordinates. Let (f?, </>) 
be arbitrary spherical coordinates on S with 9 S [0, it] 
being the azimuthal angle and (j> € [0, 2ir] being the polar 
angle. In these coordinates, the Riemannian two-metric 
q ab on S is: 

q = qee d9®d9+ <700 d^^dcji+qBij) {d6($d(f)+d4)®d9) . (29) 

The horizon may be arbitrarily distorted; q a b does not 
have to be the standard two-sphere metric. Note that 
on a sphere, any two-form L a b can be written uniquely 
as L a b — Le a b, where L is a function on S, and e a & 
is the area two-form on S; e = ^/dctqdO A d(f> where 
det q — qoeq<t><i> — q\^ is the determinant of q a b- Any one- 
form £ a can be expanded as £ = £gd9 + ^d^. The co- 
variant derivative of a one-form £ a is expressed in terms 
of the Christoffel symbols r a h C as 



V a £& = da^b — r a f, c £ c 



(30) 



The Riemann tensor of q a b has only one independent 
component 



Rabcd — -^R^ab^cd ■ 



(31) 



Now we must choose a closed curve in order to find the 
Killing vector at a single point. The equator (9 = n/2) 
is a convenient choice for the curve since it avoids the 
coordinate singularity at the poles. The tangent vector 
v a is then simply d$. The equations (|28|l then become 



(In these formulas we do not sum over repeated indices.) 
These are three coupled, linear, first-order differential 
equations in (£g, L). The same equation holds for any 
line of latitude (9 = constant). The second-order Runge- 
Kutta method was used to solve these equations. The 
initial data required for this equation are the values of 
(?0, at, say, = 0. The solution of the equation 
will be (£e, £</,,£) at cj) = 2ir. We are eventually inter- 
ested only in (£#,£0), but the function L is necessary to 
transport the data. The solution to these equations can 
be written in terms of a matrix M: 




= M 



(0=2tt) 




(33) 



(0=0) 



To find the matrix M, we start with the initial data sets 
(1, 0, 0), (0, 1, 0) and (0, 0, 1). The solutions will give the 
first, second and third columns respectively of M. Next 
we find the eigenvalues and eigenvectors of M. The eigen- 
vector with unit eigenvalue is what we want 4 . Numer- 
ically, no eigenvalue is exactly equal to unity; therefore, 
in practice, we choose the eigenvalue closest to unity (see 
the next section) . For the horizon to be axisymmetric or 
close to Kerr in any sense, this eigenvalue should be very 
close to unity. If this is not the case, then this proves that 
the horizon is not close to Kerr in any sense. All eigen- 
values will be unity in the spherically symmetric case. 

Having found the eigenvector at the point = 0, we 
then transport it to every grid point on the sphere. The 
curves used to transport the eigenvector are the lines 
of latitude and longitude. Transport along constant 9 
curves is done by equations l|32|l . while for the constant 
6 curves we use: 



86 

89 
8L 



— Tge £g + Too £,<p , 



(34) 



(Again, no summation over repeated indices.) Finally, 
having found (£g, £0, L) at each grid point, we now need 
to normalize the Killing vector £ = (£0,^) so that its 
integral curves have affine length 2ir. To do this, we 
need to follow the integral curves of £ a : 



dO 



and 



dt 



= ^{9^) (35) 



and normalize the affine parameter t so that its range 
is [0,27r]. Numerically, we only have to make sure that 



dio 



8L 



+ r e ^£0 - Ly/detq, 



(32) 



4 In principle, we should verify that every closed curve starting and 
ending at the point (8 = n/2, <f> = 0) gives the same eigenvector. 
However, as in any numerical method, we only do this for a small 
number of curves. The eigenvector obtained in this way is the 
only possible candidate for a Killing vector. 
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£ a does not vanish at the starting point. While solving 
equation i|35|) numerically, we will need the value of £ a at 
points not included in the grid. We use a second order 
interpolation method for this purpose. This finally gives 
us the normalized symmetry vector tp a , which is used to 
calculate Ja from equation (|25|l . 

IV. COMPUTATION OF J A AND TESTS OF 
THE NUMERICAL CODE 

In this section, we apply our approach to finding the 
Killing vector and our ability to identify an isolated hori- 
zon. In order to validate our approach for identifying 
Killing vectors on S, we first test our method using an- 
alytic data (jijjKij) in a simple case for which the lo- 
cation of S and its Killing vectors are also known; we 
consider the boosted Kerr-Schild solution [23 with the 
basic parameters of mass M — 1, spin a = 1/2, and a 
boost in the z-direction. The four-metric has the form 
gab = ??afc + 2H£ a £b where r\ a b is a flat Lorentzian metric, 
H is a smooth function and i a is null (see e.g. |23J for 
a detailed discussion of such metrics). The notion of a 
boost is well defined for metrics in the Kerr-Schild form 
because of the presence of the flat background metric r^. 
By boosting the black hole, we impose a coordinate dis- 
tortion on the horizon, while retaining its physical prop- 
erties. For these test cases, we know that the horizon is 
isolated; and we take advantage of only needing to com- 
pute quantities intrinsic to a two-sphere and use spherical 
coordinates. The following steps were used to test the nu- 
merical code maintaining second order accuracy at each 
step: 

1. From the Kerr-Schild data, calculate analytically 
the apparent horizon two-metric q a b, the normal 
R a , and the components of the extrinsic curvature 
K a b at the location of the apparent horizon. Dis- 
crctizc these; quantities using a. spherical grid on the 
apparent horizon. 

2. Using the discretized data, find the unnormalized 
Killing vector, £ a , at a single point (in our case we 
choose this point to be (8 = tt/2, <j> = 0)) through 
the procedure described in the previous section, ap- 
plying the Runge-Kutta method. 

3. Solve both equations O and ^ to find £ a ev- 
erywhere on the apparent horizon. 

4. Normalize the Killing vector, <fr a , using interpola- 
tion and a Runge-Kutta method for equations (|35[) . 

5. Calculate Ja via equation l|25|l . using R a given by 
the apparent horizon and <p a determined by steps 
1-4. 

The first step is easy if we have the analytic expres- 
sions for the relevant quantities. In the second step, we 
have to find the matrix M described in equation i|33[) and 




Spin parameter a 

FIG. 2: Plots of the real and imaginary parts of the eigen- 
values of the matrix M (defined in eqn. (1331 ) versus a large 
range of the spin parameter a. 



find its eigenvector with eigenvalue closest to unity. One 
can ask whether there is any ambiguity in choosing the 
right eigenvalue; is it possible that more than one eigen- 
value is close to unity? In the spherically symmetric case 
(a = 0), all eigenvalues are equal to unity and it is im- 
material which one we choose; the angular momentum 
will be zero. When a is sufficiently large, one eigenvalue 
is much closer to unity in magnitude as compared to the 
other two. In our case, it turns out that the matrix M 
has one real eigenvalue A and two complex eigenvalues 
A_R e ± iXi m . Figure [3 a plot of the real and imaginary 
parts of the eigenvalues as a function of a (for the un- 
boosted Kerr-Schild hole) , demonstrates the unambigu- 
ous nature of the eigenvalue for large values of a. The 
ambiguity may arise when a is very small. In figurc|21 we 
plot both functions for a smaller range of a. Both plots 
were generated for a resolution of dip = ir/80. The fig- 
ures show that the correct eigenvalue is typically easy to 
identify, because the other eigenvalues diverge from unity 
rather rapidly and also, at least in this case, the 'wrong' 
eigenvalues are complex while the correct eigenvalue is 
real. 

Having found the correct eigenvector and therefore the 
Killing vector at a single point, we then find it at every 
other grid point and use it to calculate Ja . Figure 0] 
plots the values of the angular momentum of the black 
hole found using equation (|25|l versus different values of 
the boost parameter for the Kerr-Schild data. Three 
different resolutions are plotted, showing a second-order 
convergence rate towards the known analytical value of 
Ja = 0.5 as expected. Although there is a slight loss 
in accuracy as the boost approaches the speed of light, 
the angular momentum loses only 1% in accuracy for 
the least resolved case (dip = 7r/20) in figure 0] when the 
boost parameter is increased from to 0.8. This loss in 
accuracy is purely due to numerical effects and is smaller 
for the higher resolution cases. We also obtained similar 
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Spin parameter 1 a 
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FIG. 3: Plots of the real and imaginary parts of the eigenval- 
ues of the matrix M versus the spin parameter a, when a is 
small. 



d 4> = 71/80 



d <j> = TC/40 



Boost parameter V z 

FIG. 4: The numerically computed angular momentum of 
the black hole at different boosts for a black hole with mass 
M = 1 and spin a — 1/2. Three different resolutions d(j> are 
shown. 



results for boosts in other directions. 

A more realistic situation is to compute Ja and Ma 
during a numerical simulation of a black-hole spacetime, 
in which the spacetime data will be given on a spatial 
grid. To test our method in this case, we again use 
boosted Kerr-Schild data, but this time we start with 
numerical data discretized on a Cartesian mesh on a spa- 
tial slice; this mesh will not coincide with the spherical 
mesh on the apparent horizon. We use an apparent hori- 
zon finder to locate the apparent horizon S and its unit 
spacelike normal R a , and construct a spherical grid on 
the apparent horizon. Let dx and g?</> be the grid spac- 
ing of the Cartesian and spherical grid respectively. We 
want the two grids to be of similar spacing, i.e. we choose 
d(j) such that d<f> ~ dx/ R where R is the coordinate ra- 
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FIG. 5: Resolution tests for the angular momentum Ja of the 
horizon. The scenarios II - IV are explained in table [I] 



dius of the apparent horizon. The data are then inter- 
polated onto the spherical grid, an additional source of 
error. We extract the two-metric q a i, numerically from 
the data, use it to find the Killing vector field (p a , and 
apply our formula for angular momentum. We present 
a series of test cases involving the black hole in boosted 
Kerr-Schild data. One is static, and three others have 
a spin of 1/2 about the z-axis, with one of the spinning 
holes boosted perpendicular to the spin along the x-axis, 
another parallel to the spin along the z-axis. Table|Ulists 
the different scenarios. Due to the additional complexity 
of having the data in a mesh that is not the one on S 
where the calculation is done, we have to deal with two 
different numerical grids. We refined both the Carte- 
sian grid and the spherical grid intrinsic to the apparent 
horizon to perform convergence tests. All runs were per- 
formed with three resolutions: 1. dx — 1/4, d<\> = 10°; 2. 
dx = 1/8, d<j) = 5°; 3. dx = 1/16, d<p = 2.5°. All lengths 
are measured in units of mass M where M is the mass 
of the black hole; Ja is measured in units of M 2 . Fig- 
ure [5] shows Ja versus resolution, and figure displays 
Ma versus resolution, showing second-order convergence 
to the known solutions for each of the cases described in 
table [I] In addition to Ja and Ma, we also monitor 
how well we converge to a truly isolated horizon, one in 
which the shear a is zero. Figure \7\ plots the value of a 
versus resolution and demonstrates second-order conver- 
gence toward zero. As expected due to additional errors 
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FIG. 6: Resolution test for the mass Ma of the horizon. The 
scenarios II - IV are explained in table Q 
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FIG. 7: This graph shows the L2 norm of (J(t)- We see that it 
converges to zero, indicating that the horizon is isolated. The 
scenarios II - IV are explained in table Q 



in the apparent horizon tracker, the convergence factors 
are not as good as in the case of analytic data; but are 
still acceptable for second order convergence. 



V. DISCUSSION 

In this section we want to compare our method of find- 
ing the mass and angular momentum of a black hole in a 
numerical simulation with other methods that are com- 
monly used. 

Note that the method proposed in this paper has three 
advantages: (i) it is not tied to a particular geometry (like 
the Kerr geometry), (ii) it is completely coordinate inde- 
pendent, and (iii) it only requires data that is intrinsic to 
the apparent horizon. The commonly used alternatives 
for calculating mass and angular momentum do not share 



all three of these features. 

The reader might object to (i) on the grounds that 
our formula for Ma (equation is also motivated 

by the corresponding Kerr formula. In this regard, the 
following remarks are worth mentioning. Even though 
our formula for Ma is based on the corresponding Kerr 
equation, much less use is made of the specific properties 
of a Kerr black hole than in most conventional meth- 
ods such as the great circle method described below. In 
particular, while deriving equation (|27|l . we do not make 
any assumption regarding the near horizon geometry. All 
that is required is that the members of the Kerr family 
are elements of the phase space under consideration (see 
[III Il2| for a detailed specification of the phase space) . 
The Kerr solutions are used only to provide a reference 
point for values of the black hole energy; this is somewhat 
analogous to choosing the zero of a potential in classical 
particle mechanics. In principle, we could also have cho- 
sen a different two parameter family of solutions but we 
choose the Kerr family due to its well known importance 
in black hole theory. 

Owing to the uniqueness theorems of classical general 
relativity, it is commonly believed that a black hole that 
has been created in a violent event will radiate away all its 
higher multipole moments and settle down to form a Kerr 
black hole near the horizon. One strategy for assigning a 
mass and an angular momentum to a black hole is then 
to identify the member of the Kerr family one is dealing 
with and to read off the corresponding mass and angular 
momentum parameters. 

While this strategy is physically well motivated, and 
one does expect the final black hole to be close to Kerr 
in some sense, we can refine this strategy considerably. 
The main difficulty with this method is that there are 
many subtleties and open questions regarding the issue 
of uniqueness of the final black hole. To briefly illustrate 
this point, let us consider the isolated horizon describing 
the final black hole. The intrinsic geometry of the hori- 
zon is, by definition, time independent. However, it is not 
necessary that four dimensional quantities evaluated at 
the horizon must also be time independent. For instance, 
using the Einstein equations at the isolated horizon, it 
turns out that the expansion and shear of the inward 
pointing normal n a (see equation |T[l) may not be time 
independent; these quantities decay exponentially |l3j . 
This means that the four-geometry in the vicinity of a 
NEH is generically not time independent, and hence may 
not be isometric to a Kerr solution. It is also not clear 
whether the four geometry tends to the Kerr geometry 
as we approach future time like infinity. Clearly, we need 
to prove a statement which says: 'sufficiently far' into 
the future, the spacetime metric in a 'sufficiently small' 
neighborhood of the horizon, is 'close to the Kerr metric'. 
In this statement, each phrase within quotes is ill defined. 
To prove or disprove such statements, it is imperative, at 
the very least, that one does not start with the assump- 
tion that the final black hole is Kerr. In particular, it is 
desirable that we not use geometrical properties of the 



13 



Kerr horizon to calculate angular momentum and mass. 

If one nevertheless makes the assumption that the final 
black hole is described by a Kerr geometry, one has to 
find a way to identify the particular member in the Kerr 
family. The method that is most commonly used is the 
great circle method which is based on properties of the 
Kerr horizon found by Smarr . It can be described as 
follows: 

In the usual coordinates, let L e be the proper length 
of the equator of the Kerr apparent horizon and L p the 
proper length of a polar meridian. Here the equator is the 
great circle of maximum length and a polar meridian is a 
great circle of minimum length. The distortion parameter 
6 is then defined to be [L e — L p ) / L e . Smarr then showed 
that the knowledge of S, together with one other quantity 
like area, L e , or L p , is sufficient to find the parameters 
m and a of the Kerr geometry. 

The difficulty with this method, apart from relying 
overly on properties of the Kerr spacetime, is that no- 
tions such as great circles, equator or polar meridian etc. 
are all highly coordinate dependent. If we represent the 
familiar two-metric on the Kerr horizon in different co- 
ordinates, the great circles in one coordinate system will 
not agree with great circles in the other system. The two 
coordinate systems will therefore give different answers 
for M and a as calculated by this method. In certain 
specific situations where one has a good intuition about 
the coordinate system being used and the physical sit- 
uation being modelled, this method might be useful as 
a quick way of calculating angular momentum, but it is 
inadequate as a general method. 

The problem of coordinate dependence can be dealt 
with in axisymmetric situations without assuming that 
the coordinate system used in the numerical code is 
adapted to the axial symmetry. The idea is to use the 
orbits of the Killing vector as analogs of the lines of lati- 
tude on a metric two-sphere. The analog of the equator is 
then the orbit of the Killing vector which has maximum 
proper length. This defines L e in an invariant way. The 
north and south poles are the points where the Killing 
vector vanishes, and the analog of L p is the length of a 
geodesic joining these two points (because of axial sym- 
metry, all geodesies joining the poles will have the same 
length). Since this geodesic is necessarily perpendicular 
to the Killing vector, we just need to find the length of a 
curve which joins the north and south poles and is every- 
where perpendicular to the Killing orbits. With L e and 
L p defined in this coordinate invariant way, we can follow 
the same procedure as in the great circle method to cal- 
culate the mass and angular momentum. This method 
can be called the generalized great circle method. 

How does the generalized great circle method compare 
to our method? From a purely practical point of view, 
note that this method requires us to find the Killing vec- 
tor, to determine the orbit of the Killing vector with 
maximum length, and to calculate the length of a curve 
joining the poles which is orthogonal to the Killing or- 
bits. The first step is the same as in the isolated horizon 



method presented in this paper. While the next step 
in the isolated horizon method is simply to integrate a 
component of the extrinsic curvature on the horizon, this 
method requires more work, and furthermore, the numer- 
ical errors involved are at least as high as in the isolated 
horizon method. Thus the simplicity of the great circle 
method is lost when we try to make it coordinate invari- 
ant, and it retains the disadvantage of relying heavily on 
the properties of the Kerr geometry. 

It should also be mentioned here that there exist ex- 
act solutions to Einstein's equations representing static, 
non-rotating and axisymmetric black holes. Examples of 
such solutions are the distorted black hole solutions found 
by Geroch and Hartle or the solutions representing 
black holes immersed in electromagnetic fields |25j . The 
apparent horizons in all these solutions are distorted, and 
the generalized great circle method will give a non-zero 
value for the angular momentum. Our method based 
on equations (|25|) and (|22|) will give Ja = because the 
imaginary part of \&2 vanishes at the horizon for these so- 
lutions. While these solutions are probably not relevant 
for numerical simulations of binary black hole collisions, 
they represent physically interesting situations in which 
a black hole is surrounded by different kinds of external 
matter fields which distort the black hole; these black 
holes may have some relevance astrophysically. These 
solutions show that the generalized great circle method 
cannot be correct in general. They also illustrate that the 
great circle method will in general give results that are 
different from the ones obtained using Komar integrals. 

A completely different approach to finding the mass 
and angular momentum of a black hole in a numerical 
solution is to use the concept of a Killing horizon. Since 
in a numerical simulation one is interested in highly dy- 
namical situations, one can not assume the existence of 
Killing vectors in the whole spacetime. Instead one as- 
sumes that stationary and axial Killing vectors exist in a 
neighborhood of the horizon, and then uses appropriate 
Komar integrals to find the mass and angular momen- 
tum. 

While this method is coordinate independent and 
does not rely on a specific metric, it has two disadvan- 
tages when compared with the isolated horizon approach. 
First, it is not a priori clear how the stationary Killing 
vector is to be normalized if it is only known in a neigh- 
borhood of the horizon. Secondly, this method requires 
the Killing vectors to be known in a whole neighborhood 
of the horizon. Computationally this is more expensive 
than finding a Killing vector just on the horizon. Concep- 
tually it is also unclear how big this neighborhood of the 
horizon should be. Furthermore, at present there is no 
Hamiltonian framework available in which the boundary 
condition involves the existence of Killing vectors in a fi- 
nite neighborhood of the horizon. In a sense, the isolated 
horizon framework extracts just the minimum amount of 
information from a Killing horizon in order to carry out 
the Hamiltonian analysis and define conserved quantities. 
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VI. APPLICATIONS AND FUTURE 
DIRECTIONS 

One situation where the calculation of Ma and Ja 
might be useful is in studying properties of initial data. 
Consider, for example, an initial data representing a bi- 
nary black hole system. If the black holes are far apart, 
then we may consider them to be isolated (this can be 
verified, for example, by calculating the shear of the in- 
ward null normal as in equation (|20[l h Then the differ- 
ence E = M^ + M A — Madm represents the binding- 
energy between the holes plus the energy due to radi- 
ation and possibly also the kinetic energy of the black 
holes. In order to find the individual black hole masses 
M A and M A equations J2HI and fTTf will be useful for 
this purpose. It is also interesting to compare different 
initial data sets which represent roughly the same physi- 
cal situation by calculating the quantity E. A lot of work 
in this direction has recently been carried out by Pfciffcr 
et al. ptlj . The isolated horizon framework may provide 
some further insights. 

The isolated horizon framework may also be used to 
construct initial data representing two (or more) black 
holes far away from each other. We want to specify the 
individual black hole spins, velocities and masses in the 
initial data when the two black holes are very far apart. 
For this purpose, the formulae for Ja and Ma would 
be relevant, and we may assume that the black holes 
will be isolated at least for a short time and will form 
an isolated horizon. The isolated horizon conditions will 
then yield boundary conditions at the apparent horizon 
which can be used to solve the constraint equations on 
the spatial slice. In fact, pioneering work in this direction 
has already been carried out by Cook [27|: the quasi- 
equilibrium boundary conditions developed by Cook are 
identical to the isolated horizon boundary conditions in 
many ways. 

Let us also briefly discuss one more important future 
application: extracting radiation waveforms. We ex- 
pect that information about gravitational radiation will 
be encoded in a component of the Weyl tensor such as 
^4 = C a bcdn a rn b n c rn d where n a and m a are members of 



a null tetrad. In certain situations where we have a back- 
ground metric or if we have Killing vectors, there may be 
a natural choice of the null tetrad used to calculate ^4. 
However, in general situations, it is not clear how to con- 
struct this preferred null tetrad. Without a well defined 
way of calculating radiation waveforms, it is very hard to 
even compare the results of different simulations which 
model similar physical situations but in different coordi- 
nates. It turns out that the isolated horizon framework 
could be used to provide a preferred coordinate system 
in the vicinity of an isolated horizon. For this purpose, 
we need more structure on the horizon than provided by 
a NEH; we need the notion of an isolated horizon as dis- 
cussed briefly at the end of section ITTaI It can be shown 
that, generically, an isolated horizon has a preferred fo- 
liation which may be said to define its rest frame. In 
general, this preferred foliation will not agree with the 
foliation given by the apparent horizons in a particular 
choice of Cauchy surfaces in spacetime. However, this 
preferred foliation can be constructed in a completely 
coordinate independent manner. Given this preferred fo- 
liation, we can construct a coordinate system analogous 
to the Bondi coordinates constructed near null infinity 
|lOl Il3| . This preferred coordinate provides an invariant 
way of calculating radiation waveforms and comparing 
the results of different simulations. 
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